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Within a mean field approximation, we study a nontopological soliton solution of the Polyakov 
quark-meson model in the presence of a fermionic vacuum term with two flavors at finite temperature 
and density. The profile of the effective potential exhibits a stable soliton solution below a critical 
temperature T < for both the crossover and the first-order phase transitions, and these solutions 
are calculated here with appropriate boundary conditions. However, it is found that only if T < Tf[, 
the energy of the soliton Mm is less than the energy of the three free constituent quarks 3Mq. As 
T > T^, there is an instant delocalization phase transition from hadron matter to quark matter. 
The phase diagram together with the location of a critical end point (CEP) has been obtained in T 
and pL plane. We notice that two critical temperatures always satisfy < T^. Finally, we present 
and compare the result of thermodynamic pressure at zero chemical potential with lattice data. 

PACS numbers: 12.39.Fe,12.39.Ba,12.38.Aw,11.10.Wx 


I. INTRODUCTION 

It is widely believed that at sufficiently high temperatures and densities there is a quantum chromodynamics (QCD) 
phase transition between normal nuclear matter and quark-gluon plasma (QGP), where quarks and gluons are no 
longer confined in hadrons[lj0. The study of the QCD phase transition is experimentally supported by the heavy-ion 
collisions at ultrarelativistic energies, accomplished in the most sophisticated accelerators, such as the Relativistic 
Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory and the Large Hadron Collider (LHC) at CERN. 
These conducted experiments provide us with the opportunity to infer fundamental information about which type 
phase of matter, hadronic or quark-gluon plasma, is stabilized in the various regimes. In order to explore a wider range 
of the QCD phase transition up to several times the normal nuclear matter density, the new Facility for Antiproton 
and Ion Research (FAIR) at Darmstadt, the Nuclotron-based Ion Collider Facility (NICA) at the Joint Institute for 
Nuclear Research (JINR) in Dubna and the Japan Proton Accelerator Research Complex (J-PARC) at JAEA and 
KEK, will make such extreme conditions possible through collisions Q. Therefore, it will be possible to test the 
theoretical predictions about the hadron-quark phase transition at high density but moderate temperature. 

On the theoretical side, the property of confinement which becomes relevant at large distances or equivalently low 
energies, has hindered the development of analytical and numerical methods capable of describing the low-energy 
nonperturbative cases, especially if baryons are involved. Therefore, the challenge for nuclear physicists remains 
to find models that can bridge the gap between the fundamental theory and our wealth of knowledge about low 
energy phenomenology. Moreover, these models should be successful in explaining empirical facts at low energies, 
for example, the dynamical breaking of chiral symmetry and the confinement, which are both intimately related 
to the nonperturbative structure of the C^D vacuum. To mention a feWjthese effective models are the MIT bag 
modelj^, the Nambu-Jona-Lasinio (NJL)j5[@ and the linear sigma modelQ for quark matter, while others are the 
Brueckner-Hartree-Fock (BHF) theory^ and the relativistic mean-field (RMF) modelsQ for nuclear matter. 

Since the strongly interacting matter at very high energy should have quarks and gluons as the degrees of freedom, 
while nucleons and mesons are the relevant degrees of freedom in the hadron phase. On one hand, the phenomenology 
of the hadron-quark phase transition is often studied with the above mentioned effective models (lol - [l^ . where the 
BHF theory including the realistic baryon baryon interaction or the RMF models are used to describe hadron phase, 
while on the other, the quark phase is treated as a thermodynamic bag model or simulated by the NJL model (or 
its modernized version, the Polyakov Nambu-Jona-Lasinio model) (l^ . The motivation for these studies, aims at 
rendering two models with different degrees of freedom compatible. In more detail, the MIT bag model and PNJL 
models do describe the very well known properties of quark matter, but they fail to reproduce the bulk nuclear matter 
and finite nuclei properties. On the other side, the BHF theory and the RMF models which are constructed for 
nuclear matter are often questionable when extending to investigate high density regimes as it is common for neutron 
stars. Then the approach of considering the mixed phase of hadron and quark matter becomes important, and the 
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proper equation of state (EOS) for the hadron-quark phase transition could be derived based on the Gibbs conditions 
for phase equilibrium. Unfortunately, this kind of study is applicable only if the QCD phase transition is of first 
order. According to a recent study on QCD phase diagram based on chiral effective models, including quark and 
meson fluctuations via the functional renormalization groupfl^, the biggest part of the QCD phase diagram shows 
crossover transitions rather true first-order ones. Thus, it is necessary to search for an alternative effective model which 
can provide a proper description of hadron-quark phase transition beyond the first-order transition cases, including 
the crossover situations, while keeping at the same time the correct degrees of freedom in quark phase and hadron 
phase intact. The nontopological soliton model rooted on the Polyakov Quark Meson Model[^ appears to fulfill the 
requirement. 

The Polyakov Quark Meson Model has so far facilitated the investigation of the full QCD thermodynamics and phase 
structure at zero and finite quark chemical potential, while it has been shown that the bulk thermodynamic predictions 
of the model agree well with the Lattice QCD data fl^ - l^ . On the other hand, starting from the same Lagran gian , 
bound states (solitons) of valence quarks can be constructed through the interaction with a and tt mesons 2J, l25l |. 
Such a nontopological soliton give rise to nucleons in the hadron phase. Moreover, the nontopological soliton model 
has been proven to be a successful approach for the description of the static properties of nucleons in vacuumjH, [23 - 
[ 3 ^. Combining these two features together, while also requiring a soliton embedded in a thermal medium, the model 
provides a suitable working scheme to simultaneously study both the restoration of chiral symmetry and the possible 
dissolution of the soliton, which simulates the deconfinement transition of nuclear matter to quark matter. 

In fact, the nucleon has been previously investigated in Ref[3ll| by employing the chiral soliton model and viewing 
it as a U = 1 chiral soliton in a cold quark medium. However, the parameters /^, and Wo- were chosen to be 
the medium-modified meson values within the NJL model. For finite temperature, Abu-Shady and Mansour have 
studied nucleon properties [s^l by employing the one-loop phenomenological mesonic potentialand the coherent- 


pair approximation[29|[30|. Furthermore, the nucleon properties as well as the thermodynamics of the system both 
at finite temperature and density are examined in Refs [3J-[^. However, these studies based on the chiral soliton 
model or other nontopological soliton models [s^ - l^ suffer from two problems: the one is that they only predict a 
first-order phase transition and the other is that the critical temperature is extremely low (Tc ~ llOMeV) as compared 
with lattice data. In this work, we will improve theses previous studies by combing the chiral soliton model with 
the Polyakov loop field. Such an extension will allow us to inspect both the crossover and first-order QCD phase 
transitions and compare directly with the lattice QCD simulations. 

The structure of the paper is as follows: in the next section we introduce the PQM model with two quark flavors. 
In Sect. HI, after obtaining the effective potential in the mean field approximation, we explore the possible stable 
soliton solutions in the model. Section IV is devoted to derive the equations of motion of the nontopological soliton 
model both in vacuum and at finite temperature and density. Section V contains the static properties of nucleon at 
finite temperature and density and the phase diagram at T —^ plane. The study of the hadron-quark phase transition 
is presented in section VI. We conclude with a summary and discussions in Sec.VH. 


II. THE MODEL 


We will be working in a generalized Lagrangian of the quark-meson model for Nf = 2 quarks and W = 3 color 
degrees with quarks coupled to a spatially-constant time-dependent background gauge field remesenting Polyakov 
loop dynamics (the Polyakov-quark-meson model or the PQM in short). The Lagrangian reads |la| 


C = tjj 


g{(T + Z75T • 7?)] V’ + 2 {dfj.ad'^a + cI^tt • 


- U{a, ir) -W($,$*,r). 


( 1 ) 


Here, we have introduced a flavor-blind Yukawa interaction of strength g, coupling the isodoublet spin-^ quark fields 
ijj = ( m , d), with the spin-0 isosinglet a and the isotriplet pion field if = ( tti , 772 , 773 ). In addition, there exists a spatially 
homogeneous time-dependent gauge field represented by the Polyakov loop potential. The coupling of the quarks with 
the uniform temporal background gauge field is implemented through the covariant derivative — lA^ and 

the spatial components of the gauge fields have vanishing background = (i^oAo- 
The purely mesonic potential for the a and tt is defined as 


U (a, if) 



-Ha- 



( 2 ) 


and the minimum energy occurs for chiral fields cr and if restricted to the chiral circle in the physical vacuum: 


■TT =/^, 


( 3 ) 
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where /,r = 93MeV, corresponds to the pion decay constant and = 138MeV is the pion mass. The last two 
constant terms in Eq.^ are used to guarantee that the energy of the vacuum in the absence of quarks is zero. The 
constant H is fixed by the PCAC relation which gives H = f-Krn^. 

The quantity is the Polyakov-loop effective potential. The Polyakov loop field $ is defined as the 

thermal expectation value of the color-trace of the Wilson loop along the temporal direction 

$ = (TreL)/iVe, = (TiyLt)/jv^. (4) 


The Polyakov loop L is a matrix in color space and explicitly given by 


L{x) = Vexp 


i / drA4(x, r) , 
Jo 


( 5 ) 


with j3 = l/T being the inverse of temperature and A 4 = iA^. In the so-called Polyakov gauge, the Polyakov-loop 
matrix can be given as a diagonal representation [4l|. Within this diagonal representation, $ and $* are complex 
scalar fields. Their mean values are related to the free energy of a static, infinitely heavy test quark (anti-quark) at 
spatial position x. The Polyakov loop expectation value (<i>) vanishes in the confined phase where the free energy of a 
single heavy quark diverges, while in the deconfined phase it takes a finite value since the center symmetry becomes 
spontaneously broken [4^ . 

The temperature dependent effective potential ,T) is constructed to reproduce the thermodynamical be¬ 

havior of the Polyakov loop for the pure gauge case in accordance with lattice QCD data, and it has the Z{3) center 
symmetry like the pure gauge QCD Lagrangian. In the absence of quarks, we have $ = <!>* and the Polyakov loop is 
taken as an order parameter for deconfinement. For low temperatures, U has a single minimum at $ = 0, while at 
high temperatures it develops a second one which turns into the absolute minimum above a critical temperature Tq, 
and the Z{3) center symmetry is spontaneously broken. The simplest Z{3) symmetric polynomial form based on a 
Ginzburg-Landau ansatz is proposed in Ref.(43| 

( 6 ) 

with 

h{T) = ao + ai + a2 ' + ag . (7) 

A precise fit of the constants Oj, bi is performed to reproduce the lattice data for the pure gauge theory thermodynamics, 
as also the behavior of the Polyakov loop as a function of temperature. The corresponding parameters are 


oo = 6.75, ai = -1.95, 02 = 2.625, 

03 = -7.44, 63 = 0.75, 64 = 7.5. 


( 8 ) 


Originally, the critical temperature Tg for deconfinement in the pure gauge sector is fixed at 270 MeV, in agreement 
with the lattice results. However, in fully dynamical QCD, fermionic contributions and the matter backreaction modify 
the pure gauge potential to an effective glue potential, which carries a flavor and chemical potential dependence of 
Tq. The actual value of Tg for two quark flavors is Tg = 208MeVfTB| (Tsj . 

A convenient framework of studying phase transitions is the thermal field theory. Within this framework, the finite 
temperature effective potential is an important and useful theoretical tool. In this section, in order to investigate the 
temperature and the chemical potential dependence of the nontopological soliton, let us consider a spatially uniform 
system in thermodynamical equilibrium at temperature T and quark chemical potential p,. In general, the grand 
partition function reads 


Z 


Trexp[-(H - M7)/r] 

/ / 'Dljj'D'ipexp 


{C -b 


( 9 ) 


where J^ = i dt fy d^x, V is the volume of the system and ^ = /r_B/3 for the homogeneous background field. 

We evaluate the partition function in the mean-field approximation similar to the work of [43| . Thus we replace 
the meson fields by their expectation values in the action. In other words, we neglect both quantum and thermal 
fluctuations of the meson fields. The quarks and antiquarks are retained as quantum fields. The integration over the 
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fermions yields a determinant which can be calculated by standard methods [1^. This generates an effective potential 
for the mesons. Finally, we obtain the thermodynamical potential density as 

r!(r,/r) = ^^ + (10) 

with the quarks and antiquarks contribution 

Siw = I jgfE, - 2N,T J 0; [l„s,+ + l„s-] , (11) 


where, Nf = 2, Nc = 3 and Eg = is the valence quark and antiquark energy for u and d quarks, and 

the constituent quark (antiquark) mass Mq is defined as Mq = ga^ together with ay = The first term of 

Eg. (HID denotes the fermion vacuum one-loop contribution, regularized by the ultraviolet cutoff. In the second term 
g^ and g~ are dehned as taking trace over color space 


4 

9q 


\ _l_ 3(4* -\- X e _|_ g 3(i?g 

\ _|_ 3((j)* -(- ^ g^-(Eq+E}/T _|_ ^-d,{Eq+iJ,)/T 


( 12 ) 

(13) 


The fermion vacuum one-loop contribution is frequently omitted. In this work, we shall include the effect 
of vacuum fluctuation on the thermodynamics. This term can be properly renormalized by using the dimensional 
regularization scheme as done for the two-flavor case in Refs. fisl - f^ . and the renormalized contribution of the fermion 
vacuum loop reads 


ll--, = = - 


NyNf 

87r2 




(14) 


where A is the arbitrary renormalization scale. It is worth to note that the thermodynamic potential and all physical 
observable are independent of the choice of A, and the A dependence can be neatly cancelled out by redefining the 
parameters in the model. 

Now the first term in the right hand side of Eq. dm, describing the vacuum contribution, will be replaced by the 
appropriately renormalized fermion vacuum contribution as given in Ea. (ll4D . Accordingly, the thermodynamic grand 
potential in the presence of appropriately renormalized fermionic vacuum contribution in the Polyakov quark meson 
model will be written as 


llMF(T,/i,CT„$,$*) = 11 m(ct«)+W($,$*,T) + X> (15) 

Here, for convenience we define a new mesonic potential 

11m(ct„) = t/(cr,7f)(16) 

which is independent of the temperature T and the chemical potential g. Minimizing the thermodynamical potential 
in Eg. lfTSl) with respective to ay, $ and $*, we obtain a set of equations of motion 

cAImf 9Hmf 9Hmf 

The set of equations can be solved for the helds as functions of temperature T and chemical potential g, and the 
solutions of these coupled equations determine the behavior of the chiral order parameter dy and the Polyakov loop 
expectation values $, as a function of T and g. 

There are two values of the constants left in the model which we need to fix, namely, and g. Unlike the pion 
meson, the mass of the sigma meson still has a poorly known value, but the most recent result of the Particle Data 
Group considers that rua can vary from 400 MeV to 550 MeV with full width 400 — 700 MeV[4^. The coupling 
constant g is usually fixed by the constituent quark mass in vacuum within the range of 300 ^ 500 MeV, which gives 
g ~ 3.3 ~ 5.3. In this work we take = 472 MeV and g = 4.5 as the typical values. It has been proved in Ref.js^ 
that this set of parameters can describe the properties of nucleon in vacuum successfully. 

In figure [U the temperature dependence of the normalized chiral order parameter dy and the Polyakov loop ex¬ 
pectation values $, 3>* at ^ = OMeV and g = 320MeV are shown in relative units. The temperature behaviors of 
the chiral condensate and Polyakov loop condensate show that the system experiences a smooth crossover transition 
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FIG. 1: (Color online) The normalized chiral order parameter cr„ and the Polyakov loop expectation values as functions 

of temperature for /r = OMeV and /r = 320MeV. The solid curves are for /r = OMeV and the dashed curves are for /r = 320MeV. 

at zero chemical potential, while there is a first-order phase transition for larger chemical potential because both the 
chiral order parameter and the Polyakov loop expectation values make jumps across the gap of the condensates near 
the critical temperature. Traditionally, the temperature derivative of the chiral condensate for u and d quarks has 
a peak at some specific temperature, which is established as the critical temperature for the chiral phase transition 
for both the crossover and the hrst-order transitions. Thus, for zero chemical potential, the chiral restoration occurs 
at ~ 201MeV, whereas for a relatively larger chemical potential ^ = 320MeV, the critical temperature moves to 
the lower temperature region around ~ lOSMeV. 

Although different from the chiral phase transition, we are still not in a position to conclusively identify the 
deconfinement phase transition through the Polyakov loop expectation values $, <1>* or their temperature derivatives 
The temperature derivatives and $*' do show one peak or more peaks for zero chemical potential or finite 
chemical potentials in calculations, but unfortunately these peaks are fake signals for defining the critical temperature 
of the deconfinement. In the next section, based on the effective potential at finite temperature and finite chemical 
potential, we will explain that there is no obvious clue to define the critical temperature of the deconfinement phase 
transition simply by using the Polyakov loop expectation values $, or their temperature derivatives if T < Tq, 
even in the first-order transition region. This is a serious problem which appeared already in the Polyakov quark 
meson model[l^[2l| or the Polyakov-Nambu-Jona-Lasinio model[4l|[43, and still persists in more recent theories. In 
the following, based on the nontopological soliton model, we will provide a distinct clarification on this point which 
will allow us to propose a convincing definition of the deconfinement critical temperature. 


III. EFFECTIVE POTENTIAL AND NONTOPOLOGICAL SOLITON 

The basic ideas behind the nontopological soliton are best illustrated by considering the original model: the 
Friedberg-Lee modelji^ or its descendant modelsjl^, [H, In these models, the confinement of quarks is ap¬ 
proximated through their interaction with the phenomenological scalar field, ct, which is introduced to describe the 
complicated nonperturbative features of the QCD vacuum. In mean field approximation, the a field has a bag like 
soliton solution, named as nontopological soliton. This in contrast to topological defects which are stabilized by the 
topological properties of the vacuum manifold. The existence of this kind of solution is closely related to a potential 
describing the nonlinear self-interactions of the a field. In general, the potential leading to the soliton solution has 
three extrema: one local minimum corresponding to a perturbative vacuum state located at cr ~ 0 , one absolute 
minimum corresponding to a physical vacuum at its vacuum value (f„, and a local maximum lying between 0 and 
dy. Therefore, the soliton solution has a spherical cavity-like structure: at large radius r, the a field assumes its 
vacuum value ay, but at small r, the a field has a value close to the second minimum of the potential near zero. In 
the Friedberg-Lee model the quarks interact with a mean a field only, this means that, in the physical vacuum state 
the quark mass is more than IGeV which makes it energetically unfavorable for the quark to exist freely, so that the 
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effective heavy quarks have to be confined in hadron bags. Sometimes it is also called as “absolute” confinement, 
similar to the MIT bag model. However, as it is known for the chiral soliton model, pions can produce strongly 
attractive forces among the quarks. Including mean pion fields also allows the meson fields to remain close to the 
minimum of the mexican hat potential, and so quarks possess physical constituent masses in the physical vacuum. 
The state is to be bound if only the total energy of system is lower than the energy of three free constituent quarks 
in the system, making thus transparent for considering the chiral soliton as a bound state in this work. 

Normally, the self-interaction potential of the u field is chosen to have a quartic form in the nontopological soliton 
model, and the coefficients in the quartic potential can be chosen so that they belong to the three typical forms as 
described as in a seminal work done by Goldflam and Wiletsfs^. In their work, they have shown that in order to ensure 
the stability of the two vacuum states and guarantee the existence of the stable soliton solution, it is indispensable for 
the potential of the a field to exhibit three distinct extrema. In the following discussion, this will be also considered 
as a key criterion for determining whether there exist stable soliton solutions for the mesonic fields and the Polyakov 
loop fields. Then with employing the thermodynamic grand potential in the presence of appropriately renormalized 
fermionic vacuum contribution in Ea. (jl5|l . we can explore the possible nontopological soliton solutions owned by the 
model under this criterion. 

However, in contrast to the chiral soliton and the Friedberg-Lee models, in the present study there are three order 
parameter variables (7„, <I> and d>* in the grand canonical potential Hmf in Eg. (1151) . so it is extremely difficult to 
investigate and demonstrate the effective potential at finite temperature and chemical potential via evolving these 
variables simultaneously in such a large parameter space. In order to simplify the problem and provide a more intuitive 
insight into the physics, we separate the study into two cases: (1) the mesonic field direction, and (2) the Polyakov 
loop field direction. For the first case, we treat the a field as a variable in the grand canonical potential Hmf while 
fixing the Polyakov loop fields on their expectation values $ nd i>*. In contrast, for the second case, the Polyakov 
loop d> and d>* are consider as variables while the a field will maintain its expectation value all the time. 




o (MeV) a (MeV) 


FIG. 2: (Color online) (a) The grand canonical potentials Hmf as a function of the chiral order parameter a for p = OMeV by 
fixing the Polyakov loop on their expectation values. (b)The grand canonical potentials Hmf as a function of the chiral order 
parameter a for p = 380MeV by fixing the Polyakov loop on their expectation values. Hmf is scaled by a factor of 

The first case is shown in Figl^l where the left panel is for zero chemical potential and the right panel for p = 
380MeV. Here, the expectation value of the pion field is chosen in the standard way as (tt) = 0. For p = OMeV, one 
clearly observes a smooth crossover of the symmetry breaking pattern. The energy difference Ae between the global 
minimum and the local maximum of the potential decreases upon the increase of the temperatures. When a critical 
temperature ~ 20IMeV is reached, Ae vanishes, which indicates that the chiral symmetry is restored. Moreover, 
according to the above criterion for the existence of the stable soliton solution, for zero chemical potential, we can 
find the stable soliton solutions at various temperature from zero temperature until to the critical temperature for 
the chiral phase transition T^. The result is believed to be held for all crossover transition region in the QCD phase 
diagram. 

For [I = 380MeV, one clearly observes the characteristic pattern of a first order phase transition: two minima 
corresponding to phases of restored and broken symmetry are separated by a potential barrier and they will become 
degenerate at T = T“. Chiral symmetry is approximately restored for T > T°, where the minimum at perturbative 
vacuum cr ~ 0 becomes the absolute minimum as shown in the right panel in Fig|^ The bag constant B is now 
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negative, then it is physically prohibited to support the existence of the stable soliton solution, so that there is no 
soliton solution anymore. Therefore, we can only obtain the stable soliton solution for T < T°, and this is also 
applicable for the whole first-order transition region in the QCD phase diagram. Moreover, the barrier between the 
two local minima of the effective potential around T“, shown in the right panel in Figl^J will decrease with decreasing 
/i. At a specific chemical potential, the barrier will finally disappear and the transition will become of second 
order. The point C (T°, of the phase diagram will be termed as the critical end point (CEP). 




O (or O’) O 

FIG. 3: (Color online) (a) The grand canonical potential IImf as a function of the Polyakov loop O (or O*) for p = OMeV 
by fixing the chiral order parameters on their expectation values. (b)The grand canonical potential Hmf as a function of the 
Polyakov loop O for /i = 380MeV by fixing the chiral order parameter on their expectation values. IImf is scaled by a factor 
of r’‘. 

Let us now investigate how the grand canonical potentials IImf evolve with the Polyakov loop field $ for different 
chemical potential, by fixing the chiral order parameters on their expectation values. The scaled grand canonical 
potential is shown in Figl3]as a function of $. From Figl3]it is obvious that these effective potentials share similar 
behaviors for both fjL = OMeV and fjL = 380MeV: there is only one minimum for each of the effective potentials, and 
these minima correspond to the expectation values $ of the Polyakov loop field at specific temperature and density. 
In the crossover transition region, with the raising of the temperature, the expectation value <l> moves to its higher 
value smoothly and slowly. But in the first-order transition region at high density, accompanied by the jump of chiral 
order parameter a, the expectation value $ develops a disconnection across the gap from relatively small value to its 
maximum, which indicates that there exists a degenerate value of the Polyakov loop variable $ along the first-order 
transition line in the QCD phase diagram for very high chemical potential. This implies that the integration of the 
quark and meson fields in the grand canonical potential Hmf would only result in a trivial Polyakov loop effective 
potential at finite T and /r, and such a naive potential does not tell anything about the critical point at which the 
deconfinement transition should definitely happened. Therefore, the jump induced by the chiral order parameter is 
not to be supported by the effective potential of the Polyakov loop field itself, then it certainly can not be treated as a 
signal for the deconfinement phase transition. This is the reason why we argue that there exists no obvious criterion for 
defining the critical temperature for the deconfinement phase transition in terms of using the Polyakov loop variables 
$, <i>* or their temperature derivatives, as long as the temperature T is smaller than the critical temperature Tq for 
deconfinement in the pure gauge sector. 

Nevertheless, the advantage is that we do not have a bag-like soliton solutions for the Polyakov loop variables $, 
$*, since there is only one minimum in the effective Polyakov loop potential. On the other hand, the Polyakov loop 
variables $, will always develop their expectation values $ and in whole space, such that these fields should 
be regarded as homogeneous background thermal fields on top of which the chiral soliton is to be added on. 


IV. NONTOPOLOGICAL SOLITON SOLUTION IN THE MODEL 

In vacuum, the Polyakov loop variables $, 3>* are set to zero and the thermodynamic grand potential Hmf reduces 
to the purely mesonic potential Dm. Following reference [s^, in the mean field approximation, the a and tt are taken 
as time-independent, classical c-number fields, which only differ from their vacuum values in the neighborhood of the 
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FIG. 4: (Color online) The quark fields in relative unit and the a, tv fields scaled with /,r as function of the radius r in vacuum. 


quark sources. The state of the quarks {^„(r)} with energy {e„} and the cr(r), 7r(r) meson helds satisfy the coupled 
set of the Euler-Lagrange equations of motion 

-ia-V0 „(r)-g^[(T(r)+f75r-7r(r)](/)„(r) = e„(/)„(r), (18) 

-V^g(r) + ^„(r)(/)„(r) (19) 

^occ 

-V^7r(r) + =-g^ ^„(r)i 75 T()i„(r) (20) 

^occ 

with 

j 4i{^)<Pn{v)(f’r =l ( 21 ) 

where a and /3 are the conventional Dirac matrices. 

The ground state of the chiral soliton is the state with N quarks in the same lowest Dirac state (pQ, with energy 
e. In the following, our discussions are constrained in the case of = 3 for baryons. In order to obtain solutions of 
minimum energy, we adopt the “hedgehog” ansatz with the meson fields are spherically symmetric and valence quarks 
are in the lowest s-wave level 


a = a{r), tt = r7r(r), 


4‘o 


<r) \ 

ia ■ rv(r) I 


where x is a state in which the spin and isospin of the quark couple to zero: 


( 22 ) 

(23) 


(<? + r)x = 0. 


(24) 


Now the system is spherical symmetric and the Euler-Lagrange equations of motion (fT8l) - (|20l) transform in radial 
coordinates to 


du{r) 

dr 


dv(r) 

dr 

d^a{r) 

dr'^ 


2 da{r) 

4 - 7 — 

r dr 


d?‘TT{r) 2 dirir) 


dr"^ 


-k 


dr 


■■-ie + ga{r)) v{r) - gn{r)u{r), 

(25) 

- 57r(r)^ v{r) -k (e - gcr{r)) u{r), 

(26) 

-Ng{u^{r) v'^{r)) , 

(27) 

27r(r) dVlM ^ ^ ^ 

^2 g^ = ^Ngu{r)v{r), 

(28) 
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and the quark functions should satisfy the normalization condition 


47 r J + u^(r)) dr = 1. 


These equations are subject to the boundary conditions which follow from the requirement of finite energy: 

^(0) = 0,^^ = 0,7r(0) = 0 
dr 

u{oo) = 0, cr(oo) = /^, 7r(oo) = 0. 


(29) 


(30) 

(31) 


The asymptotic vacuum value of the soliton field has to be determined by an additional condition, i.e. that the 
physical vacuum is recovered at infinity. In this “physical” vacuum the quarks are free Dirac particles with the 
constituent mass gay , and the chiral symmetry is spontaneously broken. By solving the coupled differential equations 
(l25l) . (l26l) . (lT7l) and (IMI) with the normalization and the appropriate boundary conditions, in Fig|4]we plot the a, tt 
and quark helds prohles in arbitrary unit as functions of r for zero temperature and chemical potential. 

If we put N quarks into the lowest state with energy e, the total energy of the hedgehog baryon is given by 


E = Ne + Att 


11 
2 UE) ^2 


diT 

dr 


+ 


+ ^^m(o'i,) 


dr, 


(32) 


which is normally identified as the mass of the nucleon below. 

As a next step we consider a B = 1 localized bound state (soliton) in a thermal medium. Customarily, the thermal 
medium can be treated as a quark medium or a nuclear medium due to the interaction of the three valence quarks 
with the quark Dirac sea and the nucleon Fermi sea via the meson fields [5ll|. 

For the quark medium, the thermal medium is filled with quarks of a constituent mass Mq, and the soliton energy 
is given by the sum of the energy of the valence quarks, the meson fields and their interactions as shown in Eg. (1321) . 
Then a new set of coupled equations of motion for the meson helds could be derived by simply replacing the relevant 
mesonic potential IIm with the thermodynamic grand potential Hmf- Accordingly, a set of coupled equations for 
mesons can be described as 


d^air) 2 da{r] 

dr"^ r dr 


diluF 


d‘^'K(r 

dr'^ 


2 dirir) 

A - -1 — 

r dr 


da 

27r(r) cIOmf 


= Ng (u^(r) - v'^ir)) , 
= —2Ngu{r)v{r). 


dn 


(33) 

(34) 


For satisfying the requirement of Hnite energy of the soliton, one of the boundary conditions in Eq. (I5T|) should modihed 
accordingly as: r —>■ oo, a{r) approaches to the expectation value dy, where thermodynamic grand potential Hmf has 
an absolute minimum. 

As long as the unbound constituent quarks, treated as the homogeneous background thermal helds with T and 
/i, are allowed to penetrate into the soliton, they will bring an additional contribution to the total baryon density. 
Thus to ensure that the solitonic baryon number is equal to one, the normalization condition equation (1291) should be 
modihed as 


47r J r^ {u^{r) + v^{r)) dr = 1 — Bm, 


with 


Bm = 47r 


/ P^r^dr. 
Jv 


(35) 


(36) 


Here, and V being the volume of the soliton. 

In contrast, when a soliton is embedded in the medium of nucleons we have to consider a Fermi sea of nucleons 
instead of quarks, due to conhnement. This is because the Dirac sea consists of quarks and therefore only determines 
the vacuum sector. Thus in this case the mesons are directly coupled to the nucleons of the Fermi sea. Accordingly, 
the terms representing the thermal medium effects in Eas. (l33B4p should be modihed as: 


dfl 


MF 


da 

dfiuF 

dn 


gNi'fpN^’N), 

gN{4’Nij5T1pN)- 


(37) 

(38) 
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Here, the bracket () denotes the expectation value of the operator between the nuclear ground state, ipN is the nucleon 
field and gijv is a coupling constant which relates the nucleon mass to the non-zero expectation value of the scalar 
meson field. Unluckily, both the coupling constant qm and the scalar and pseudoscalar densities of nucleons and 
antinucleons in the above equations cannot be obtained from the present model self-consistently. In fact, they have 
to be considered as input parameters which could be taken from the RMF approximation or the BHF theory 00- 
Therefore, in order to avoid inconsistencies, it is customary to treat a hot and dense thermal medium as a uniform 
constituent quark and gluon medium (quark medium) with solitons embedded, as in (m|(51|[52|. 

As long as there is no bag-like soliton solution for the Polyakov loop variables in whole space, the Polyakov 

loop variables $, take constantly their expectation values $ and $*. Hence these variables denote the contributions 
only to the thermodynamic grand potential Umf rather than to the equations of motion for the nontopological soliton 
solutions. Consequently, the properties of a soliton placed in a thermal medium can be investigated by solving the 
four coupled Euler-Lagrange equations that arise from the thermodynamic grand potential Umf in Ea. dT5|) . This 
system of equations does not possess analytic solutions, but is readily solved numerically. Various numerical packages 
are available for the solution of such equations. One that has been widely used in this field, is COLSYSjl^. 


V. NUCLEON STATIC PROPERTIES AT FINITE TEMPERATURE AND DENSITY 

We firstly study soliton solutions at finite temperature and density by solving the coupled differential equations 
(ESI) . (E51) . (1551) and (IMl) with the normalization condition and the appropriate boundary conditions. In FiglSl we plot 
the u{r), w(r), u{r) and v{r) fields at zero and finite chemical potential = 380 MeV) for different temperatures. 
These two chemical potentials correspond to the typical crossover and first-order phase transitions in the QCD phase 
diagram, respectively. For both cases, it is shown that all the fields are moving towards the trivial values while the 
temperature increasing. When T is lager than some critical temperature T^, there only exist trivial solutions for 
the coupled equations of motion and solitons are melted away. These trivial solutions indicate the restoration of the 
chiral symmetry in full space. Moreover, the lack of solitonic solutions are usually considered as a signal for the 
delocalization of the baryonic phase. 




r{fm) 


r(fm) 


FIG. 5: (Color online) (a) The quark fields in relative units and the a, n fields scaled with f-^ as function of the radius r for 
T = 177MeV as ^ — OMeV. (b) The quark fields in relative units and the cr, tt fields scaled with as function of the radius r 
for T = 56.1MeV while /r = 380MeV. 


Based on the above analysis, for both the crossover and first-order transitions, the effective potential supports the 
existence of the stable soliton solution for the meson fields, as long as T is lower than T^. This implies that the 
baryonic phase can be indeed found in the chiral symmetry-breaking phase. However, the stability of such baryonic 
phase should be checked carefully by comparing the total energy of the system in thermal medium with the energy of 
three free constituent quarks in the system. By subtracting the homogeneous medium contribution. the total 
energy of system Mn is plotted as a function of the temperature for = OMeV and ^ = 380MeV in FiglHl Here, 
one finds that for a smooth crossover of the symmetry breaking patten at low density region, both Mn and 3Mq fall 
smoothly from the corresponding vacuum value as T goes to high temperature. When T is close to some high critical 
temperature simeq2QVbAeN for fi = OMeV, both Mjv and 3Mg experience a steep descent region. But as shown 
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FIG. 6: (Color online) The total energy of system Mjv and the energy of 3 free constituent quark 3Mq as function of the 
temperature T. Here one set is for /r = OMeV, another set is for /r = 380MeV. 


in Figl in the high temperature regime 3Mg drops more quickly than that of as T > 177MeV, 3Mq < Mm- 
This implies that even though the stable soliton solution still exists in the temperature region T € [177, 201]MeV, it 
is energetically unfavorable. The baryonic phase will definitely melt away into three free constituent quarks. In this 
way, we can identify this specific temperature ~ 177MeV as a critical temperature of the deconfinement phase 
transition. 

Based on the right panel in FiglHl one can see that in the large density region the total energy Mm decreases 
monotonically with increasing the temperature T from zero to higher values. As the temperature approaches the 
critical temperature T(:, Mm starts to deviate from the ones in vacuum signihcantly. When T > T^, Mm jumps to 
zero quickly, which indicates the delocalization phase transition from nucleon matter to quark matter due to the fact 
that the effective potential does not support the existence of the stable soliton solution. The energy of three free 
constituent quark 3Mg (or ay) shows the similar behavior as Mm- By comparing the two energies in Figl®! we can 
show that for T < the nucleon bound sate is stable and 3Eq is larger than Mm, but the difference decreases with 
the increase of temperature, and the two energies begin to cross over at the critical temperature T°. Therefore, the 
critical temperature for the deconfinement phase transition is coincident with that of the chiral phase transition for 
the first-order phase transition. 

We infer the occurrence of the chiral phase transitions of u and d quarks and the deconfinement phase transition 
at finite temperature and finite density, and show the T — fj, phase structure of the Polyakov quark meson model in 
FigE] based on the nontopological soliton picture. For two light flavors, there is a crossover in the low density region 
and a first-order phase transition in the high-density region, and in the middle position there exists a critical end 
point (CEP). From figure [71 the critical temperature for the deconfinement phase transition is lower than that 
of the chiral phase transition, and both two critical temperatures decrease smoothly as n goes to high value. With 
the increasing of /i, the difference between and becomes smaller and smaller, while at some critical chemical 
potential it becomes zero, which identifies the critical end point (CEP) for a second order phase transition. The 
corresponding values are (T°,/r^) ~ (119, 302)MeV. 

Here are several remarks on the phase diagram presented in PigjTl Prom the above discussions on the effective 
potential and nontopological soliton concerning the deconfinement and chiral phase transition, we conclude that the 
effective potential always support the existence of the stable soliton solution in the system if T < r°, but not for 
T > T°. It is required that the critical temperature defined as the deconfinement phase transition in a nontopological 
soliton model usually less than the critical temperature for the chiral phase transition, as < T^. In this study, in 
the first-order region, we take the “=”. In contrast, in the crossover region we should have the “<”. This conclusion 
is in qualitative agreement with the result shown in Pig.6 in Ref.[l^ at relatively low and middle densities. But for 
high density, they produce a very strange behavior for the deconfinement crossover phase transition for the Polyakov 
loop variables $, <i>*. This is different from the Eriedberg-Lee model [37l - [^ and its descendant modelsfs^liSl, which 
only predict a first-order phase transition in the phase diagram. The chiral soliton model combined with the Polyakov 
loop has an obvious advantage in the description of QCD phase diagram, since it can allow the prediction of the 
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H (MeV) 


FIG. 7: (Color online) Two-flavor phase diagram in the T — fj, plane in the Polyakov quark meson model based on the 
nontopological soliton pictnre. The dash-dotted curve is the critical line for which characterizes the confinement phase 
transition, and the dashed lines are the critical line for conventional chiral phase transition in the region of crossover. The solid 
lines indicates the first-order phase transitions, and the solid circle indicates the critical end points for chiral phase transitions 
of u and d quarks. 



FIG. 8: (Golor online) The proton charge ‘r.m.s’ radius of a stable chiral soliton as a function of temperature T at p = 0 MeV 
and fj, = 380 MeV. The solid curve is for fj, — 0 MeV while the dash-dotted curve for p = 380 MeV. The dotted curve is for the 
unstable baryonic phase existing in the crossover phase transition when > T > . 


crossover transition at low and middle chemical potential. Nonetheless, it shows the first-order phase transition for 
high chemicalpotential. The result is in agreement with other predictions demonstrated in effective models and lattice 
QCD data[ll-I^. 

At the end of this section, the proton charge r.m.s radii i? of a stable chiral soliton as a function of temperature for 
/i = 0 MeV and ^ = 380 MeV are illustrated in FiglSl it gives a signal of a swelling of the nucleon when temperature 
and density increase. In both cases, R increases slightly at low temperatures while the latter is increased. As T 
approaches T^, R sharply grows and disappears. Another interesting result displayed in FiglS] is that the maximal 
radius R at various densities are almost same when T near T^, which hints that solitons start to overlap each other 
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with the similar expansion rate at for different densities. 


VI. QCD THERMODYNAMICS AT ZERO CHEMICAL POTENTIAL 


In order to investigate the influence of the Polyakov loop on the equilibrium thermodynamics of the system, we 
calculate the pressure of the system P during the QCD phase transition form hadron phase to quark phase according 
to two different models as follows. 

Firstly, we adopt the mean-field approximation as usual by replacing a, tt and the Polyakov loop variables $, 
<i>* by their expectations values. In other words, we neglect both quantum and thermal fluctuations of the meson 
fields and the Polyakov loop variables but retain the quarks and antiquarks as only quantum fields in the entire phase 
diagram. This is of course not a realistic scheme, especially at low T and /r, since due to the confining forces quark and 
antiquarks will recombine into mesons, baryons and antibaryons. Hence, the character of the chiral phase transition 
described by the mean-field approximation could be drastically changed in hadronization process from quark phase to 
hadron phase. But, if we discard these affects, all thermodynamic quantities can be obtained from the grand canonical 
potential in a spatially uniform system Hmf in Eq- P5I) , which is determined as the logarithm of the partition function. 
The negative of grand potential which is normalized to vanish at T = ^ = 0 gives the thermodynamic pressure in 
PQM model 


PpQM — 


(39) 


The pressure obtained in the above equation can be directly compared with lattice data. 

However, the hadron and quark phases can be distinguished by empirical facts and phenomena at low and high 
energies. At low temperature and low baryon density, the hadronic phase exhibits a dynamical breaking of chiral 
symmetry and the confinement, and the baryon and meson act as the active degrees of freedom here. On the contrary, 
at very high temperature or baryon density, quarks and gluons will be set free to play the dominant roles in quark 
gluon plasma. Such a scenario can be realized in the nontopological soliton (NS) model vividly as follows: in the 
hadron phase, the state of the free quarks is not the ground state of the strongly interacting matter, and as a result 
three valence quarks will form the bound state of the nucleon. Therefore, the hadron phase only possesses baryons and 
mesons. On the other side, when T > T^, the solitons are going to dissolve, and the hadronic phase will eventually 
evolve to quark phase with free quarks. 

For simplicity, within the NS model, we assume an ideal case of the system by taking the hadronic phase as 
a noninteracting hadron gas composed of nucleons and tt, a mesons with the effective masses M^, and 

in thermal medium. The Polyakov-loop variables are treated as the background thermal fields and accordingly the 
Polyakov-loop potential has been subtracted already. It is then straightforward to write down the normalized pressure 
of the system in terms of nucleons and mesons for the hadronic phase [^(45j| 


pH 

^NS 





(40) 


Where vn = 4 for nucleon, z/,r = 3 for pion and t'o- = 1 for sigma meson. The last term B*{Mn) is introduced in 
order to recover the thermodynamical consistency of the system, since the nucleons are treated as the chiral solitons 
with a temperature-dependent masses [s^. The explicit expression of this term can be evaluated by the additional 
constraint {BPep/ dMM)T = 0, which gives 


B* {Mn{T)) 


= B* (Mjv(O)) - z/jv 



dT' 


MNiT') 


d^p 1 


1 


1 


,(41) 


with E'jy = \Jj? + Mm{T')^. The energies in Eq.(j40|) Em = \J = \JfP + M.n.(T)^ and E^ = 

\JfP' + are corresponding to nucleon, pion and sigma mesons, respectively. Mm is obtained as the energy of 

soliton, whereas, the a and tt masses are determined by the curvature of Omf in Eg. (HSU at the global minimum: 




(42) 


Since it is unfavorable for solitonic nucleons to survive at high energy when the temperature is across the deconfine¬ 
ment critical temperature ^ 177MeV, the baryonic bound state formed by three constituent quarks will definitely 
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FIG. 9: (Color online) The normalized pressures variation with respect to temperature for PQM model and NS model ai fj, — 0 
with ~ 201MeV. The solid line corresponds to the pressure in PQM model, while the dash-dotted line is for the pressure in 
nontopological soliton model. All calculations are compared to lattice data {Nt = 6) from Ref. [s^. 


dissolve and the system should now be regarded as a quark phase including the free quarks, mesons and the gluons 
mimicked by the Polyakov loop. Consequently, the pressure of the NS model in terms of free quarks and mesons 
incorporating with the Polyakov-loop potential in quark phase results in 


pQ _ j. 


T J ^ {in [l + + In [l + } 


-v^T 


[ d^P f 

J (27r)3 I 


In 


1 — e 


-E^jT 


} - ".r I {l„ [l - I _ J,). 


(43) 


Here, Vq = 2NcNf = 12 is the number of internal degrees of freedom of the quarks and Eg = yjp + is the valence 

quark and antiquark energy for u and d quarks, and the constituent quark (antiquark) mass Mq is given by Mq = gSy. 

Thermodynamic pressures divided by the QCD Stefan-Boltzmann (SB) limit are illustrated at /r = OMeV in Figl9] 
for three models, and the pressure has already been normalized to vanish at T = /i = 0. For Nf massless quarks and 
Nl — 1 massless gluons in the deconfined phase, the QCD pressure in the SB limit is given by 






— (—Y -\- - . 

180 127r2'r 


& 


(44) 


where the hrst term is the gluonic contribution and the rest involves the fermions. 

At very high temperature, e.g. around twice the chiral critical temperature, the pressure of the PQM tends to 
approach that of the NS model in the quark phase. This implies that there only exists a weak interaction between 
particles in quark phase, and the quasi-particles model is a good approximation for the description of the weak- 
interaction quark gluon plasma. However, when the temperature decreases, the Ppqm deviates from the pressure 
more and more. When the T reaches out to the chiral critical temperature ^ 201MeV, the gap between the two 
pressures arrives at its maximum value and then decrease smoothly until zero temperature is reached. For comparing 
with lattice simulations with a temporal extent of Nr = 6 which is closer to the continuum limit in Ref. , the strong 
interaction between the particles in the NS model cannot be simply discarded and it really plays a dominant role in 
producing a correct and reasonable thermodynamical pressure of the system. 

In Fig El we also plot the pressure as a function of the temperature starting from the hadron to the quark phase at 
/is = 0 for the NS model, while varying the baryon masses for various temperature and densities in the confined phase. 
From the figure, the dash-dotted curve shows rapidly changed discontinuities at the deconfinement critical temperature 
from hadron matter to quark matter. This indicates a pseudo first-order phase transition for the delocalization 
transition and signals a drastic structural change for nucleon when the system goes with the diffusions of the solitons 
(nucleons) into thermal medium simultaneously. This strange behavior of the pressure Pns around the deconfinement 
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critical temperature is believed to be removed by bringing in a self-consistent interaction of cr, w, p mesons with 
nuclear matter (solitons) in hadron phase. The strong interactions between mesons and nucleons, which are widely 
adopted in nuclear matter and finite nuclei would provide necessary suppressions on the pressures for the 

mesons and nucleons in hadron phase, and the remaining self-interaction of the mesons in quark phase could give 
further suppression on the pressure around T“, similar to the PQM model. Finally, it is worth to note that, besides 
the interactions, another important influence which was neglected here is the center-of mass (c.m.) correction to the 
nucleon properties, which have been considered in soliton bag models and the QMC model largelv [6ll - [^ . These two 
corrections are certainly out of the scope of our current topic and we prefer leaving them for a future study. 

Unlike the case in the quark phase, here the pressures of the PQM model and the nontopological soliton model are 
suppressed in the hadron phase and start to increase when deconfinement sets in. The small difference among the 
pressures of the Pns and Fpqm is due to the different treatment of the mesons in the hadron phase. In the NS model 
the mesons are taken as the active degrees of freedom, but for the PQM model they are purely mean fields and should 
be restrained to the except values when T is close to the T^. From FiglHl it is shown that the effective nucleon mass 
Mm slightly deviates from its vacuum value with temperature increasing, only at the critical temperature T^, Mm 
experiences a sharp jump. Consequently, the contribution of nucleons to the total pressure in Ea. (|40)) is very small in 
hadron phase as far as the chemical potential is small. To estimate, it only gives 8.6% contribution to Pns when T is 
around ~ 177 MeV for zero chemical potential. On the contrary, if three bound quarks are set to be free in quark 
phase, the valence quarks will give a dominant contribution to the pressure in Eg. (1431) . 


VII. SUMMARY AND DISCUSSION 

In the present paper we have investigated possible nontopological soliton solutions in the effective potential of the 
PQM model in the presence of renormalized fermionic vacuum. The results show that as long as the temperature 
is not larger than the chiral critical temperature T^, there exist truly stable soliton solutions in the model for both 
crossover and first-order phase transitions. 

Even though there are stable soliton solutions for the Euler-Lagrange equations of the model at finite temperature 
and density, the stability of the solitons (nucleons) have to be checked and analyzed carefully in thermal medium by 
comparing the effective masses of nucleons with the energies of three free constituent quarks. Our results show that 
the chiral phase transition and the delocalization phase transition from nucleon matter to quark matter take place 
simultaneously for the first-order phase transition. For T < the free constituent quarks are not the ground state 
of the strongly interacting system, and the quarks will reorganize so to form lower-energy bound states carrying the 
hedgehog configuration. However, as soon as the temperature T crosses over the such bound states cannot survive 
anymore, and the system experiences a first-order hadron-quark phase transition to the chirally symmetric phase. 

The situation differs from the case of the first-order transition, in the crossover transition, even though the effective 
potential genuinely ensures the stable soliton solution in the system, but it is energetically unfavorable for nucleons to 
exist when T is across the deconfinement critical temperature T^. The difference between and is about 24MeV 
for zero chemical potential, but it will decrease to zero as p. increases to some value around 302MeV. This particular 
point sometimes is denoted as the critical end point (CEP), which appears as well in the phase diagram of Figl?]. 

In order to compare our results with the lattice QCD simulations and other models at zero chemical potential but 
finite temperature directly, we have investigated the thermodynamic properties of the nontopological soliton model in 
the PQM model. It is found that the inclusion of the Polyakov loop is necessary and important when comparing with 
the lattice QCD simulations. When compared with the previous studies in this field, we notice a quite improvement 
on the topic by providing a reasonable critical deconfinement temperature ~ 177MeV for p = 0 and extracting a 
standard QCD phase diagram in agreement with the lattice data and other phenomenological models’ predictions [3|. 
However, the strange behavior of the performed pressure in the hadron-quark phase transition indicates that the 
description of the hadron phase as a non-interacting hadron gas of the nucleons and mesons with medium-modified 
masses has underestimated the important effects of their interactions, and these interactions should be introduced 
to further suppress the pressures of the mesons and nucleons both in hadron and quark phase. In other words, the 
present form of study is a prototype and still not suitable for the proper description of nuclear matter and finite nuclei 
in hadron phase yet, and it deserves further efforts on making the model applicable for hadron-quark phase transition 
completely and satisfactorily. Eventually, work in this direction is in progress. 
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